Effect of interference on thermal noise and coating optimization in dielectric mirrors 
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Optical multilayer coatings of high-reflective mirrors significantly determine the properties of 
Fabry-Perot resonators. Thermal (Brownian) noise in these coatings produce excess phase noise 
which can seriously degrade the sensitivity of high-precision measurements with these cavities, in 
particular in laser gravitational- wave antennas (for example project LIGO), where at the current 
stage it is one of the main limiting factors. We present a method to calculate this effect accurately 
and analyze different strategies to diminish it by optimizing the coating. 

Traditionally this noise is calculated as if the beam is reflected from the surface of the mirror 
fluctuating due to the sums of the fluctuations of each layer. However the beam in fact penetrates 
a coating and Brownian expansion of the layers leads to dephasing of interference in the coating 
and consequently to additional change in reflected phase. Fluctuations in the thickness of a layer 
change the strain in the medium and hence due to photoelastic effect change the refractive index of 
this layer. This additional effect should be also considered. It is possible to make the noise smaller 
preserving the reflectivity by changing the total number of layers and thicknesses of high and low 
refractive ones. We show how this optimized coating may be constructed analytically rather then 
numerically as before. We also check the possibility to use internal resonant layers and optimized 
cap layer to decrease the thermal noise. 

PACS numbers: 04.80.Nn, 42.79.Wc, 07.60.Ly, 05.40.Ca 
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I. INTRODUCTION 

Any precise measurement faces a challenge of different 
noises superposing useful signal. Brownian noise coming 
from chaotic thermal motion of particles is one of the 
examples. The Michelson interferometer is able to de- 
tect small changes in the lengths of its arms: two beams 
traveling different optical paths, interfere on the detec- 
tor producing intensity which depends on the difference 
between phases of the beams. Thermal (sometimes also 
called Brownian) noise in coatings and substrates of the 
mirrors of the interferometer results in fluctuations of 
mirrors' surfaces averaged over the beam spot and adds 
a random phase to the waves. This effect is one of the 
key factors limiting the sensitivity of laser gravitational 
wave detectors Though the thickness of the mul- 
tilayer coating is just several micrometers, the internal 
mechanical losses in layers is several orders of magnitude 
larger then in the substrate. That's why coating ther- 
mal noise, in accordance with the fluctuation-dissipation 
theorem, dominates and exceeds other noises produced 
in the mirrors 

In this paper we analyze different effects and strategies 
aimed to decrease the thermal coating noise for general- 
ized multilayer reflective coating. Traditionally this noise 
is calculated as if the optical beam is reflected from the 
surface of the mirror fluctuating as the incoherent sum 
of the fluctuations of each layer. However, the beam ac- 
tually penetrates the coating and Brownian expansion of 
the layers leads to dephasing of interference and conse- 
quently to additional change in reflected phase Q . Fluc- 
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tuations in the thickness of a layer change the strain in 
the medium and hence due to photoelastic effect change 
the refractive index of this layer. This additional effect 
should also be considered. It is possible to make the noise 
smaller while preserving the reflectivity by changing the 
number of layers and thicknesses of high and low refrac- 
tive components [3, H| ■ We also check the possibility of 
using internal resonant layers @ and an optimized cap 
layer to decrease the thermal noise. 

Brownian noise is not the only source of noise produced 
by the coating. Fluctuations of temperature, which 
translate into displacement of mirror's surface through 
thermal expansion (thermoelastic noise) 0,111 and change 
of optical path due to fluctuations of refraction index 
(thermorefractive noise) [9jl combine producing general- 
ized thermo-optical noise [3, G3 • Brownian fluctuations 
causing displacement of mirrors' surface and the previ- 
ously neglected correlated photoelastic effect produced 
by these fluctuations both lead to fluctuations of the op- 
tical path and should be treated simultaneously. The 
Brownian branch of noises, which is the topic of this pa- 
per, and thermo-optic one are uncorrelated as they rep- 
resent uncorrelated fluctuations of volume and tempera- 
ture. 



II. PHASE NOISE FROM MULTILAYER 
COATINGS 

A. Reflectivity 

To calculate the amplitude and phase of a reflected 
beam the impedance method [11] will be used below. We 
found this method more convenient for analytical consid- 
eration than the equivalent and more widely used matrix 
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method [lj. 

To consider reflection on each boundary separating the 
layers starting from the substrate/coating boundary (see 
pic. [TJ we introduce an effective impedance Z[z) and 
amplitude reflection coefficient T{z) as follows: 
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where E and H are tangential electric and magnetic fields 
in the standing wave, while E + , H + and E-, H— are for- 
ward and backward (reflected) waves, n is the refraction 
index, fj, and e the relative permeability and permittivity, 
and Z v is the vacuum impedance (Z v — 1 in Gaussian 
CGS system). 

As tangential fields E and H are continuous in a 
medium without free currents, the effective impedance is 
also continuous on all boundaries, while reflection coeffi- 
cient experiences jump. Meanwhile the reflection coeffi- 
cient changes continuously between boundaries according 
to the following equation: 

Tp p iujt+ikonj{z — dj) 

T(z - d,) = , ■■ — r - TT = r{z)e~ l2k ° n > d \ (4) 

where ko = =£■ is the wave vector of the optical field in 
vacuum, A is the wavelength. This allows one to calculate 
any multilayer coating layer by layer recursively starting 
from the substrate in which the impedance is equal to the 
impedance of free substrate rj s and moving to the surface, 
turning from the reflection coefficient 1^ = T(— dj — 
0) to the effective impedance Zj = Z{— dj — 0) = 
Z(— ^ dj + 0) when facing the boundary and back after 
crossing it (see Fig. [!}. In particular it is possible to 
exclude effective impedance from calculations by deriving 
the formula for reflectance jump itself: 



where g lj = (fij = 2k njdj. 

In the case of classical reflective A/4 layers, all 
impedances and reflection coefficients are real. 



B. Interference 

We now assume that each of the layers can have varia- 
tion in thickness ddj and independent variations of its 
refraction index Srij arising from fluctuations. Both 
should be accounted in ((4|), introducing —<pj — > — tpj — 
2koSnjdj — 2konjSdj = —tpj+Aj. We also have to rewrite 
Vi 17*(1 + Si]i) m which is a consequence of re- 

fraction index change Srjj = As before, moving 

layer by layer to the surface, we expand each result into 
a series to the first order of variations Srij and Sdj. In 
this way we can build a perturbed amplitude reflection 
coefficient T' : 
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Here m is the index of the layer of interest (m = 7V+e for 
total reflectance, where ./V is the full number of layers, "e" 
in indices represents external medium and "+e" means 
Taking into account that Aj, <C 1, we can find 
the equivalent phase shift 5<p as well as modification of 
reflectivity <5r (leading to amplitude noise which cannot 
be found in traditional approach) collecting all imaginary 
and real parts noting the decomposition: Te e = T(l+e). 
Total fluctuations may arise both from layer thickness 
fluctuations Sdj (Brownian and thermoelastic noises), or 
from deviations of refraction index Srij (photoelastic and 
thermo-refractive noises). 
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FIG. 1. A schematic of a multilayer coating 



C. Photoelastic effect 

Besides surface displacement noise there should be a 
photoelastic noise in the Brownian noise branch. Pho- 
toelasticity is a phenomenon of refraction index change 
under deformation: 



Aft 



■■Piji 



(8) 



where Bi is the optical indicatrix, Uj - strain tensor, 
Pij is photoelastic tensor and indices i,j € f;6 [13| . In 
case of cylindrical symmetry we have longitudinal effect 
ABi = Pi3«3 = paSd/d and transversal effect ABi = 
Pip u pp- However, only noise produced by longitudinal 
effect has the same origin as thermal one we consider 
(movement in "z" direction), and in this way there is a 
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theoretical possibility of their interference compensation. 
For longitudinal effect of variation of refraction indices 
one obtains: 



5n : , 
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There is also nonzero Sn z component, but we consider 
only rays incoming perpendicularly to the surface and do 
not take it into account. It is known that tantalum oxide 
used in multilayer coatings Ta^O^ - is a rutile (titanium 
oxide) type crystal with tetragonal symmetry Rutile has 
P13 = 0.171, P23 = 0.16. From [3] we can make a rough 
estimation for tantalum oxide PTa 2 O s < 0.18. We put 
for simplicity P13 = P23 = PTa 2 o 5 — 0.17. The other 
component of the coating - fused silica has p i3 — p 2 3 = 
PSi0 2 = 0-27. 

The transversal effect should be considered separately 
as u pp - motion is not correlated with u zz oc Sd - mo- 
tion. This effect should not give out more noise than the 
longitudinal one and should be added incoherently 



D. Brownian branch of noises 

We use longitudinal photoelastic effect to convert fluc- 
tuations of refraction index into fluctuations thickness of 
layer: 
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where pj is effective photoelastic index for j-layer. Then 
coating induced deviations of reflected phase and refle- 
tion coefficient are 



fyc = ^2a J Sd j , 
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where 
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Let us consider one end mirror in the arm of an inter- 
ferometer. Thermal displacement of the mirror's surface 
produces phase fluctuations in the interferometer output. 
It is more intuitive to consider the case of contraction 
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FIG. 2. Phase shift of the optical wave reflecting from an 
unperturbed (upper figure) and perturbed mirror. S(fo, Sip g 
and Stp c are shift in the total phase, the shift due to the surface 
displacement and the shift due to interference dephasing in 
the coating respectively (Sdj < 0). 



(Fig. [5]) of the mirror. Then the length of additional gap 
for the light to travel before entering the mirror is — Sd 
(as Sd < for contracting), yielding the phase shift 
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The total phase shift produced by the perturbed coating 
(relative to the unperturbed one) will be 

N 

Sip* = -2k [fN + e{-l) N - j z^ jnj - 1] Sdj, (17) 

3=1 

where we took into account that inside A/4-reflector all 
quantities are real and otj = —2kQnjipjf]^ +e (—l) N ^^Zj. 

It is also important to admit that in a "good mirror" 
approximation, when 1 — |T| <gc 1 (in this case depending 
on the topmost layer Zn — > or Zn — > oo) the amplitude 
reflection coefficient correction produced by each layer 
Pj = (-l) N - 3 fN+cJjfj^ -> (for A/4-refiector). 

The term before Sdj can be treated as a noise coef- 
ficient showing a contribution of each layer into the to- 
tal noise. This coefficient can have any sign, depending 
on the values of interferential contribution ("— " sign) or 
surface displacement ("+"), but only its absolute value 
is significant as noise contributions from different layers 
are added incoherently. 

Using the acquired formulas we can plot a diagram 
of phase shift contribution of each layer and values of 
noise spectral densities in the whole. In Fig. [3] such 
distribution is plotted, keeping the sign from (|17[) . It 
can be seen that the interference part of noise plays a 
role in a few outer layers (order of penetration depth) Q 
while Brownian (surface displacement) noise forms the 
major part. Several layers can even demonstrate nearly 
complete noise compensation. 

We should note that the noise contribution of a layer 
is formally composed of three summands: main Brow- 
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43 layers 
42 layers 




FIG. 3. Noise coefficient (keeping the sign) from each layer in 
a coating consisting of 42 (circles) or 43 (squares) layers on 
silicon substrate. 



nian (surface displacement) contribution, interferential 
part and photoelastic effect: 



S(ps = 2koddj 
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drij ddj ^ ' 



(18) 



where tp denotes the phase of complex reflectivity. For- 
mulas ©-([71) give analytical expression of the derivatives. 
Their sign distribution may be illustrated as follows. If 
the coating contracts, then the phase shift produced by 
each layer is positive due to the change of its thickness 
as the Brownian (surface displacement) noise is not re- 
ally a phase shift acquired by light inside the mirror, but 
outside it (Fig. [5]). The contraction of each layer at the 
same time leads to an increase of refraction index as in 
normal materials it grows with density providing posi- 
tive phase shift. Interference dephasing (phase shift due 
to reduction of layer thickness itself), on the other hand, 
is somewhat that compensates the phase shift produced 
by both effects. It looks like photoelastic effect can play 
only negative role, however it can be the thing to reduce 
too high interference dephasing (for example it lowers the 
total noise in 42-layer A/4-coating). 

Equation (1181) is quite easy for numerical calculations 
as the partial derivatives in it may also be calculated nu- 
merically. We used this approach for independent check- 
ing of formulas ©-([7]). 



E. Noise spectral density 

Using (TTT)) one can estimate noise spectral densities if 
the noise spectral density of each layer is known. In the 
model of independent thin layers on an infinite half space 
substrate, each layer behaves just as if it was the only 
layer on the substrate. This model was heavily treated 
and the solution is well known 0, For our purpose 
however, we should split the total surface fluctuations of 



one layer into two parts. The first one being the fluc- 
tuations of the thickness of the coating layer S c and the 
second one being the fluctuations of the substrate surface 
induced by losses in the coating S s . Interference and pho- 
toelastic effects influence only the first term. If the losses 
in the layer responsible for both fluctuations (shear and 
expansion losses) are equal, what we assume in this pa- 
per, then these two spectral densities are uncorrelated. 
Otherwise cross correlation terms should be taken into 
account. This splitting may be easily obtained using 
the approach presented in [3J, assuming that the noise 
produced by each layer is independent (6d%) — > S c (fl)j, 



(SdjSdk 



0: 



(19) 

4* fl T(l + i/ i )(l-2i/ i ) 



^3 



' irw 2 n Yjil-i/j) 

Ak B TY J {l + v s ) 2 {l-2v s ) 2 
irw 2 n Y 2 {l-u 2 ) 



(20) 



where Vj is the Poisson coefficient of layer j, Yj - its 
Young's modulus (Y s and v s are the parameters of the 
substrate), <f>j is the mechanical loss angle, w is the Gaus- 
sian beam radius on the mirror, Q is frequency of anal- 
ysis, ks is the Boltzmann's constant and T is the tem- 
perature. Then we get spectral densities of phase and 
amplitude reflection fluctuations 
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The first sum can be simplified and the second is zero in 
assumption of a "good mirror" and A/4 layers: 
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S v =Ak 2 J2 [S c (ty 1 " 7 "'' 7 " 



N 



+ S s (n) m N] (23) 

for 2N layers, and S v + Ak^S\ for 2N + 1 layers, where 
a rn = for zero outer impedance (2N layers with n\ > 
n,2) and a m = n 2 n 2 for infinite outer impedance (2N + 1 
layers). 

To simplify the comparison of this type of noise with 
other types of noises and Fabry-Perot coordinate sensi- 
tivity we turn phase noise into noise of effective reflecting 
surface displacement 
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(in units of m 2 /Hz) at 100 Hz frequency. 



(24) 
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Tvne 


42 x A/4 


41 x A/4 + A/2 


43 x A/4 


Transmittance t. ppm 


2.28 


1.08 


0.54 


Brownian 10~ 20 m/\/Hz 


0.632 


0.635 


0.645 


X With interference 


1.96% 


2.34% 


1.75% 


X With photoelasticity 


2.33% 


1.85% 


1.31% 


X Modified cap 


2.33% 


2.76% 


1.81% 



TABLE I. Silica-tantala mirror efficiencies relative to the 
Brownian noise. Standard LIGO coating consists of 41 
Layers+A/2 cap mirror. Modified cap has optical width A/4 
(42 layers case). 



Calculations were made for a silica-tantala mirror of 
42-43 layers (21 pairs of S1O2 Ta-iO^ A/4-layers on fused 
silica substrate with or without additional A/4-layer). 



Vl 


= 0.17, 


ni 


= 1.45, 


Vh 


= 0.23, 


nh 


= 2.06, 


Y, 


= 7.2 x 10 10 Pa, 


<t>l 


= 0.4 x 10~ 4 , 


Y h 


= 14 x 10 10 Pa, 


<t>h 


= 2.3 x 10~ 4 



A = 1.064 x 10~ 6 m; w = 0.06 m; T = 290 K 

Resul ts w ere formed in table U in form of correction 
_ V-ggr^x/s x 100%. Numerical estimations for relative 

VbBr 

transmittance noise is St/t = 2rVSr/(l - |T| 2 ) < 10~ 12 
Hz" 1 / 2 . 

The interference correction to thermal coating thick- 
ness noise is about 6%, or 7.5% taking photoelasticity 
into account. But thickness variation of tantala layer is 
much smaller than its bending = 0.36££). That is 
why the interference correction to full coating Brownian 
(displacement) noise is only about 2.0%, or 2.3% taking 
photoelasticity into account. 

III. OPTIMIZATION STRATEGIES 

A. Additional top layer-corrector 

One may tweak the thickness of the topmost "correct- 
ing" layer in an attempt to minimize the noises using 
interference effects. This method proved to be useful for 
thermoelastic and thermorefractive noises Q ■ Using for- 
mulas ©-© we can obtain 

s v =4fc 2 £ [s c m m (t4^ - i%%T + n) 

+ s s {n) m N] + s' c 

S' c = l$t(g N +c+e)(l ±j c sin((j) c ))n c ip c - n e f S c (fl) c 
+ S S (Q) C (25) 

for 2N layers, and + 4/c§S 1 for 2N + 1 layers, where 
=K(.92Ar+c+e) and "+" for zero impedance of 




FIG. 4. Noise coefficient distribution in coating for 42 layers, 
keeping the sign. Silica substrate, vacuum medium (squares 
- simple mirror; circles - a mirror with modified layer #26 
(16 from top) d m = 0.98A/2). 

last but one layer (2iV + c layers with n-i < ri\ and a m — 

2 2 

3 * jft(ff2iv+i+c+e) an d "— " for infinite impedance of last 
but one layer (2iV+l+c layers). The index "c" represents 
the cap layer corrector and "+c" means "+1". 

Results are quite unfavorable: for even number of lay- 
ers + cap minimum of noise is at n c < 1 while its suppres- 
sion x = V s^f " x 100% is onl y °- 04% - For odd la y ers 

+ cap absolute value of noise doesn't become lower than 
6.198 x 10 - 20 m/V Hz, which means suppression is less 
than 0.69% (for n c = 3.6; d c = 0.42A/4). Even after 
removing a pair of layers, the noise is about 6.04 x 10~ 20 
m/ y/ Hz, which is more than for even number of layers. 

This means that standard coating with top silica A/2 
layer is reflectance-optimized as well as "all A/4" coating 
(cap= A/4) is noise-optimized (see Tab. 1). 

B. Layer-corrector inside the mirror. 

In J.KimbleQ proposed an idea of inserting a resonant 
layer into the mirror. This case was studied numerically 
(Fig. 2]). The maximum suppression of 4.4% was shown 
by layer-corrector close to d = A/2, which is the resonant 
cavity. But such modification increases power transmit- 
tance more than two orders of magnitude. If we add 
8 bilayers to restore transmittance, suppression will be 
more than eliminated (—14%). 

C. Two-sided and double mirror. 

A novel combined structure was propose in [lj| with 
just a few layers on the front side of a big silica substrate 
with antiresonance optical length and other layers moved 
to the the bottom (two-sided mirror or Khalili etalon). 
The idea is that only top layers produce Brownian noise, 
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FIG. 5. Noise distribution in the coating, keeping the sign. 
Silica substrate, vacuum medium. Squares -A/4 general coat- 
ing; circles - corresponding etalon with A/4 substrate; trian- 
gles - etalon, optimized for interference, (0.35A/4). 




FIG. 6. Suppression as a function of excess substrate optical 
thickness (in units of A/4) in addition to integer number of 
halfwavelengthes for different noise ratios 7^ = I 2 -. 

while bottom layers do not contribute to thermal noise 
as their surfaces move in the opposite direction. In this 
case we should also pay attention to interference effects, 
because first layers and substrate are well penetrated 
by light. That also means that coating noise and sub- 
strate noise in etalon should be treated simultaneously 
and there is a possibility of high interferential compen- 
sation (Fig. 15161) . The main difficulty with the Khalili 
etalon is its high sensitivity to the manufacturing preci- 
sion and fluctuations of its optical thickness produced for 
example by temperature variations. Namely the impreci- 
sion of substrate optical length by 0.07A/4, corresponding 
to the mirror's temperature variation of 6 mK, increases 
noise by 5%. 

The same idea may be realized in another geometry 
(double mirror or Khalili cavity) with combined end mir- 
ror consisting of two individually suspended mirrors sep- 



arated by a controlled gap. The first mirror has small 
number of layers and hence low noise, while the layers 
of the second one provide the required reflectivity. The 
sensitivity to the gap length is two times higher, though 
it may be controlled with actuators in real time yielding 
desired conditions. Our calculations show encouraging 
suppression of noise in both schemes. The deficiency of 
both schemes however is high power circulating in the 
mirror's substrate, which leads to various thermal lens- 
ing and detuning effects. 

The absolute value of maximum effect is highly depen- 
dent on ratio of thickness and bend noise spectral den- 
sities, which is yet unknown. This far we can only say 
that noise suppression and amplification effects decrease 
practically linearly with 7^ = I* (Fig. [5]). 

D. Modifying silica-tantala ratio. 

A promising way to reduce the thermal noise in the 
coating was proposed [J, [l6| , which suggests decreasing 
the thickness of lossy high-index (tantalum-oxide) layers, 
presumably preserving the total bilayer optical thickness 
to be A/2 (nidi + rihdh = A/2). To keep the required 
reflectivity, more bilayers should be used. It was found 
numerically that there is an optimum in the ratio of lay- 
ers' thickness and number of layers providing minimal 
noise at a given reflectivity. 

It appears that noise suppression x is highly dependent 
on noise ratio in layers 

Sh/dh Zh&hni . . 

x x ~Q~rr = ti — = 7 ' ( 26 ) 
ai/di imuh 

For the LIGO parameters [17] 7 = 4.56. In [4] coat- 
ing was optimized for a chosen parameter 7 = 7. An 
etalon silica-tantala mirror of 27 xA/4 layers+A/2 cap 
was numerically optimized. Resulting coating had 16 
silica-tantala bilayers with nidi = 1.3827A/4, n^dh = 
0.6173A/4, a thin cap n t d t = 0.1620A/4 and the first layer 
fihAh — 0.5560A/4 on substrate (34 layers total). In an 
experiment with this mirror design, noise suppression of 
Xexp = (8.8 ±2.0)% was observed. Our calculations with 
all material parameters taken from [4| yield Xth 7 = 8.2%, 
and if 7 = 9.23 estimated from the same experiment is 
used |li|, one gets xth = 9.1%. 

IV. OPTIMAL COATING 

It is well known that for a fixed number of bilayers, 
a multilayer coating with quarter length layers (QWL) 
with (fh = if 1 = 7r provides the largest reflectivity [l2^ . 
The LIGO interferometers however require not only large 
reflectivity but also small noise added by the coating. 
The coating consists of two different materials having 
noticeably different mechanical losses. This fact stimu- 
lated mostly numerical attempts to construct more op- 
timal coating which could have smaller noise with the 
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increased number of layers but decreased total thickness 
of the "bad" component, while still preserving the de- 
sired reflectivity [j, [l6|. It was also found empirically 
that a better coating can be achieved when tph < tt and 
it is a common knowledge that tph + ipi should be equal to 
27r. It is possible however, to construct a nearly perfectly 
optimized coating analytically and we will show that the 
"common knowledge" is in fact incorrect. 

We would like to find optimal thicknesses of the com- 
ponents of a bilayer for a given thickness ip^. Suppose 
we have a bilayer inside a coating and the amplitude re- 
flectivity on the boundary to this bilayer from the side 
of the substrate is Ti n — Toe lVo , where To is real ampli- 
tude and ipo is some initial phase. Let's introduce nota- 
tions: Ti n = IV - initial reflectivity, Ti - intermediate 
reflectivity, T out = - output reflectivity, Ti n+ i — Y%i - 
reflectivity, that will be initial for the next pair (see Fig. 
ED: 



r 2 = 



1 + 9hiT ' 
9ih + Tie-' llph 
1 + gihTxe-^h 

r 2 < = r 2 e-^<, 



(27) 

(28) 
(29) 



where gij 



rii+rij 



Now we can find the optimal phase 
Note that i^iW _ ejr^t 



ipo maximizing |r m+1| . uuav 0yo - —g^ 

does not depend on ipi. After some math we find 



tan ipo 



1 



1 



^cot^. 
+ 9 2 hl 2 

- g 2 sm((f h ). 



(30) 
(31) 



In the last approximation we used the fact that gu — 0.17 
is small. It is also important that reflectivity increases 



with new pairs of layers only if ipo £ 



(will be 



explained later). To optimize the next layer, we should 
provide the same phase <pi n for it, which gives 

Plj = fo J+1 - fo } - Phi ~ 2sin(^ 3 ) 5 L + (win), (32) 




where j stands for the bilayer number, m is integer num- 
ber. For a series of identical bilayers that means 



<Pl 



-tph - 2su\(ip h )g 2 + (irm) 



(33) 



It can be shown that in our case (ip e [— §;§], <£h € 
[0; 2tt}) m > and even. As we need to shorten layers for 
minimum noise n = 2. Now from the last equation it is 
clear that only in the case of QWL coating (pjj + ip^ = 
2tt. In other cases, however, there should be a small 
correction to maximize the reflection. Note that for the 
first layer = with indefinite phase, thus satisfying 
the requirement on ipo. 

To make analytical approach consider a bilayer some- 
where in the middle of coating. Suppose that incoming 
reflectivity is close to 1: 



ir, 



ir 



2k' 



1 



(34) 



Then, expanding formulas for |r out | into series to the 
second order of e and using (|3"Tj) we obtain 



where 



= 1 - 2ae - a(l - 2a)e 2 



(1-.9 2 ) 2 



(35) 



(gy/2(l - cob(pi)) ± ^l-2s 2 cosM+.g 4 ) 

(36) 

Here we have "+" sign, when ipo g [— ir/2; w/2]. It can be 
shown, that only in this case a < 1, meaning increasing 
reflectivity. 

Assuming \2ae\ -C 1, |(1 — 2a) e\ <C 2 we can rewrite 
local reflectivity as |r out | = 1 — ae and get total power 
transmittance in the following form: 



where 



1 + 2g e cos(ip c 



2 



■ g 2 sm{(p h )) + g 2 e 



(37) 



(38) 



describes the coating-air boundary (g e — ^~^ ! ). For 
this formula to work we need to satisfy the assumptions 
we made. Calculations for g = 0.17 give a £ [0.55; 1] 
(a(7r) = 0.55) and e <C 0.5. That also requires iph £ 
[7r/4; Tit /A]. It can be shown numerically that all those 
requirements can be satisfied with just three layers on 
the substrate. 

Now we can eliminate the total number of layers from 
equations to design an optimal mirror for a given power 
transmittance. As calculations using (J6J - (JTJ) for the 
total noise are rather complicated but provide small cor- 
rection only, we consider the simplified formula (|16[) . For 
the spectral density we obtain 



FIG. 7. A bilayer inside a multilayer coating 



Set = A((N + E)(*t<p h + pi) + tp c - (fi) 



(39) 



8 



TvDe 


25 + A/2 




Our Method 


Transmittance r, ppm 


277.5 


277.7 


277.7 


X, Brownian (displacement) 





8.16% 


8.4% 


X, With interference 


3.37% 


11.03% 


11.27% 


X, With photoclasticity 


2.63% 


10.29% 


10.5% 


Real suppression 





7.93% 


8.18% 



TABLE II. Optimization results. (7 = 7) 



in units of m 2 /Hz where A — is dimensional con- 

stant. Then for minimization we have 

— =Ejip eh , + Eip d - (f 0e + Lp + (40) 
In To — In f3 — In e 

(7^ + ipi) - tpi + ip c (41) 

In a 

Here <^ e /j f t i denote phase thicknesses of first E bilay- 
ers (from substrate), (poe, <po - initial phases for initial 
(first E layers) and regular bilayers, ip^, ipi - regular layer 
thicknesses, (f c - cap layer thickness. 

All values here may be expressed in terms of (fh, <p e h> 



L£r and minimized. The results obtained are very close to 
[J] (see table [TTJ) and represent an analytical alternative 
to numerical optimization. 

In conclusion, wc found explicit formulas for spectral 
density of phase noise produced by Brownian thermal 
fluctuations in arbitrary multilayer coatings taking into 
account interference effects and photoelasticity in the 
coating. Interference effects and photoelasticity play a 
role only in few top layers and give correction of the order 
of 2%. Some optimization methods taking into account 
inerference were considered. The method of modifying 
silica-tantala ratio was found to be the most efficient so 
far. Another promising approach is compound mirrors. 
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